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In the framework of the nonlinear A-model we investigate propagation of solitons in atomic vapors 
and Bose-Einstein condensates. We show how the complicated nonlinear interplay between fast 
solitons and slow-light solitons in the A-type media points to the possibility to create optical gates 
and, thus, to control the optical transparency of the A-type media. We provide an exact analytic 
description of decelerating, stopping and re-accelerating of slow-light solitons in atomic media in 
the nonadiabatic regime. Dynamical control over slow- light solitons is realized via a controlling 
field generated by an auxiliary laser. For a rather general time dependence of the field; we find the 
dynamics of the slow-light soliton inside the medium. We provide an analytical description for the 
nonlinear dependence of the velocity of the signal on the controlling field. If the background field is 
turned off at some moment of time, the signal stops. We find the location and shape of the spatially 
localized memory bit imprinted into the medium. We discuss physically interesting features of our 
solution, which are in a good agreement with recent experiments. 
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I. INTRODUCTION. 

Recent progress in experimental techniques for the co- 
herent control of light-matter interaction opens many op- 
portunities for interesting practical applications. The ex- 
periments are carried out on various types of materials 
such as cold sodium atoms 0, Q , rubidium atom vapors 
0OHSI! solids 0, , photonic crystals . These ex- 
periments are based on the control over the absorption 
properties of the medium and study slow light and su- 
perluminal light effects. The control can be realized in 
the regime of electromagnetically induced transparency 
(EIT), by the coherent population oscillations or other 
induced transparency techniques. The use of each dif- 
ferent materials brings specific advantages important for 
the practical realization of the effects. For instance, the 
cold atoms have negligible Doppler broadening and small 
collision rates, which increases ground-state coherence 
time. The experiments on rubidium vapors are carried 
at room temperatures and this does not require appli- 
cation of complicated cooling methods. The solids are 
obviously one of the strongest candidates for realization 
of long-living optical memory. Photonic crystals provide 
a broad range of paths to guide and manipulate the slow 
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light. The interest in the physics of light propagation 
in atomic vapors and Bose-Einstein condensates (BEC) 
is strongly motivated by the success of research on stor- 
age and retrieval of optical information in these media 

[3 ii ii ii nam. 

Even though the linear approach to describing these ef- 
fects based on the theory of electromagnetically induced 
transparency (EIT) [12| is developed in detail mod- 
ern experiments require more complete nonlinear descrip- 
tions 11]. The linear theory of EIT assumes the probe 
field to be much weaker than the controlling field. To 
allow significant changes in the initial atomic state due 
to interaction with the optical pulse, here we go beyond 
the limits of linear theory. In the adiabatic regime, when 
the fields change in time very slowly, approximate ana- 
lytical solutions 0, and self-consistent solutions 
were found and later applied in the study of processes of 
storage and retrieval [13] ■ Different EIT and self- induced 
transparency solitons of nonlinear regime were classified 
and numerically studied for their stability flal . As it 
was demonstrated by Dutton and coauthors |19| strong 
nonlinearity can result in interesting new phenomena. 
Recent experiments and numerical studies |lll l'2("| have 
shown that the adiabatic condition can be relaxed allow- 
ing for much more efficient control over the storage and 
retrieval of optical information. 

In this paper we study the interaction of light with a 
gaseous active medium whose working energy levels are 
well approximated by the A-scheme. Our theoretical 
model is a very close prototype for a gas of sodium atoms, 
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whose interaction with the light is approximated by the 
structure of levels of the A-type. The structure of lev- 
els is given in Fig. Q where two hyperfine sub-levels of 
sodium state 3 2 Si/ 2 with F — 1,F — 2 are associated 
with |2) and |1) states, correspondingly The ex- 
cited state 1 3) corresponds to the hyperfine sub- level of 
the term 3 2 P 3 / 2 with F = 2. We consider the case when 
the atoms are cooled down to microkelvin temperatures 
in order to suppress the Doppler shift and increase the 
coherence life-time for the ground levels. The atomic co- 
herence life-time in sodium atoms at a temperature of 
0.9/iK is of the order 0.9 ms Q. Typically, in the ex- 
periments the pulses have length of microseconds, which 
is much shorter than the coherence life-time and longer 
than the optical relaxation time of 16.3ns. 
The gas cell is illuminated by two circularly polarized op- 
tical beams co-propagating in the z-direction. One beam, 
denoted as channel a, is a cr~-polarized field, and the 
other, denoted as b, is a c + -polarized field. The corre- 
sponding fields are presented within the slow-light vary- 
ing amplitude and phase approximation (SVEPA) as 

E = e a £ ae i(kaz-«>at) + g b Sbe i(k b z-u b t) + c c 

Here, k a ^ are the wave numbers, while the vectors e a , e& 
describe polarizations of the fields. It is convenient to 
introduce two corresponding Rabi frequencies: 



\3) = \3P„ 1 ,F = 2,M F =0) 




FIG. 1: The A-scheme for working energy levels of sodium 
atoms. The parameters of the scheme are the following: 
lo 12 /(2tt) = 1772MHz, u/(2ir) = 5.1 • 10 14 Hz (A = 589nm), 
and A is the variable detuning from the resonance. 



that v a = v\, = vq. Hence, within the SVEPA the wave 
equations are reduced to the first order PDEs: 

dfoa = ivoPsi, d^VLb = iis p32- (4) 
Equations Eqs.Q) can be rewritten in a matrix form as 



d^Hj = i-^ [D, p] , 



(5) 



(2) In the new variables the Liouville equation takes the form 




where p a b are dipole moments of quantum transitions in 
the channels a and b. 

In the interaction picture and within the SVEPA, the 
Hamiltonian H\ = H Q + Hj describing the interaction of 
a three-level atom with the fields is defined as follows: 

Ha = -jD, H, = -i (Jl a |3)<l| + n 6 |3)<2|) + h.c, (3) 

where 



D 



Here A is the variable detuning from the resonance and 
we set h = 1. 

The dynamics of the fields is described by the Maxwell 
equations 

(0? - c 2 d 2 z )VL a e l ^ z -^ = -^d 2 t (paie 4 ^-^*') , 
(<9 t 2 - c^ 2 )^^"""*) = -737<9 t 2 (p^e 1 ^-^) , 

where v a = (n A \p a \ 2 uj a )/e , v h = (n A \p b \ 2 u b ) / e 0l n A is 
the density of atoms, and eo is the vacuum susceptibility. 
Here p is the density matrix in the interaction represen- 
tation. For typical experimental situations the coupling 
constants v a b are almost the same. Therefore we assume 



d T p : 



-D-Hj,p 



(6) 



Here £ = z/c, r = t — z/c. To make parameters dimen- 
sionless, we measure the time in units of optical pulse 
length t p = lfjts typical for the experiments on the slow- 
light phenomena gg. We also normalize the spatial co- 
ordinate to the spatial length of the pulse slowed down 



in the medium, i.e. L 



Vgtp 



Here f2o is a 



typical magnitude of the controlling field required in EIT 
experiments. We assume this field to be of order a few 
megahertz. We choose f2o = 3 as a representative value. 
This corresponds to a group velocity of several meters per 
second, depending on the density of the atoms. We take 
the group velocity to be 10 _7 c, so the pulse spatial length 
is 30/im, and £ is normalized to 10 _13 s. Then, in the di- 

mensionless units, the coupling constant vq = = 4.5. 
The retarded time r is measured in microseconds and the 
Rabi frequencies are normalized to MHz. 
The system of equations Eqs.ljSJ),© is exactly solvable 
in the framework of the inverse scattering (IS) method 
[2lL E^ . I2H . l24j . This means that the system of equa- 
tions Eqs.0!© constitutes a compatibility condition for 
a certain linear system, namely 



<9 T $ = U(X) $ = I AD $ - iHi $, 



(7) 

(8) 
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Here, A £ C is the spectral parameter. The comparison 
against leads to the zero-curvature condition 
UM U C {X) - V T (X) + [C/(A), V(X)] = 0, which holds iden- 
tically with respect to the linearly independent terms in 
A. It is straightforward to check that the resulting condi- 
tions coincide with the nonlinear equations Eqs.JHJ),©- 
At this point it is worth discussing the initial and bound- 
ary conditions underlying the physical problem in ques- 
tion. We consider a semi-infinite C > active medium 
with a pulse of light incident at the point £ = (initial 
condition). This means that the evolution is considered 
with respect to the space variable £, while the boundary 
conditions should be specified with respect to the vari- 
able r. In our case we use as the asymptotic boundary 
conditions the asymptotic values of the density matrix at 
r — > ±oo. To solve the nonlinear dynamics as described 
by equations Eqs.lJSJ),©, the IS method considers the 
scattering problem for the linear system Eq.Q, while 
the auxiliary linear system Eq.JSJ describes the evolu- 
tion of the scattering data. The purpose of this work is, 
in particular, to study an essentially nonlinear interplay 
of the fields in both the channels. This goal leads to 
considering for equation Eq . J7|) the scattering problem of 
finite density type (cf. [21j and references therein), i.e. 
Q, a ,b — > as t — ► ±oo. For an account of other re- 
sults for the A-system accessible throu gh t he IS method 
see, for example, references [bH I22I I23I l25j . In this work 
we choose to use an algebraic version of the IS method, 
i.e. the Darboux-Backlund (DB) transformations. The 
DB method does not require a full investigation of the 
initial value problem and merely allows mounting of a 
soliton on a chosen background. The resulting solution 
is, of course, consistent with the underlying initial value 
problem. We use Darboux-Backlund transformations, in 
the spirit of |2(| |27], IH, |2i| , up to certain modifications. 
In the physical case considered in this report, the system 
is assumed to be initially in a stationary state described 
by the following background solution: 

n a = o,n b = n(r), p = |V-tW«*l - |i><i|. (9) 

Notice that the state |1) is a dark-state for the control- 
ling field f2(r). This means that the atoms do not interact 
with the field f2(r) created by the auxiliary laser. The 
configuration Eq.© above corresponds to a typical ex- 
perimental setup (see e.g. 0,0,3)- The function fl(r) 
models the controlling field, which governs the dynamics 
of the system. The time dependence of this function can 
result from modulation of the intensity of the auxiliary 
laser. In general, f2(r) can al so depend on the spatial 
variable £. However, we do not specify such dependence 
explicitly in the formalism below except for a simple case 
of linear phase shift, which is discussed in section ITTll 
The paper is organized as follows. In the next section we 
describe the Darboux-Backlund transformation for the 
A-system. In section ITTT1 we describe the mechanism of a 
transparency gate for the slow-light soliton. In section llVl 
we discuss an exactly solvable example of manipulation 
of slow- light solitons, while section fVl considers a similar 



problem for the case of a fairly arbitrary controlling field. 
Section IVTl is devoted to conclusions and discussion. 



II. DARBOUX-BACKLUND 
TRANSFORMATION FOR THE A-SYSTEM 

In this section we describe the Darboux-Backlund (DB) 
transformation for the A-system. First we reformulate 
the linear system Eqs.Q,© in the matrix form, viz. 

d T ^> = l -D^ C-iH^, (10) 
g c9 = % ^L p mv. (11) 

Here $ is a matrix consisting of three linearly inde- 
pendent solutions of the linear system Eqs.J2J),ljBJ) cor- 
responding to three (not necessarily different) values of 
the spectral parameter A, i.e. A', A", A'". The matrix 
spectral parameters C is defined as 

fx' \ 
C = A" , 
\ A'" J 

while V- 1 = £ - A • I. 

The ./V-fold (N > 1) Darboux-Backlund transformation 
can be formulated as 

JV 

*[«]^(-l)» +1 V„(A)*r", Hq=/. (12) 

n=0 

It is clear that the linear system Eas. (|10|I . HllJ) is covari- 
ant with respect to this transformation provided that for 
< n < N the following Darboux-Backlund dressing 
transformations are satisfied 



Hi[N] Z N ~ n {A) = S A r_„(A)H J + id T Z N ^ n {A) 

-|[A(5jv-„+i(A)-AS iV _ ri (A))], (13) 
p[N] S N - n (A) = Hiv-n(A) P + %d ( Z N _ n+1 (A), (14) 

together with the convention Sat + i(A) = So (A) = /. 
The meaning of Eas. l|13() . l|14|l is that they connect the 
"seed" solutions H],p of the nonlinear system with the 
dressed (iV-soliton) solutions Hi[N], p[N]. To derive the 
matrices {Sfe},-,, we specify a set of solutions {^k}i-—i 
corresponding to certain fixed values of the matrix spec- 
tral parameter £, i.e. {£fc} fc=1 , where 



Afe-i 
















, P- 1 = C k - A • I 








A fe -i / 
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We then demand 



N 



(-l) n+1 S A r_„(A)* fc V- n = 0, k = 1, . . . , N. (15) 

n=0 

This linear system allows the dressing matrices 
{Sfe(A)}^ =1 to be obtained through Cramer's rule. It 
can be shown that solutions of Ea. (|15fl satisfy the rela- 
tions Eas. (P|l . (fHjl . 

Since in what follows we only discuss the case N = 1 for 
convenience, we changed notations as follows -H/[l] — > 
Hi, p[l] — > p, Si (A) — > 3(A). Then the dressing formu- 
lae Eas.l(T! ^ .!fT % lJTt |) reduce to 

Hi = Hi - \ [D, S(0)] , p = 3(A) p S _1 (A) (16) 
* = ^V- 1 -5(A)*, (17) 

while from the linear system Eq. (|15fl we obtain 

5(A) = *!(£! -A-/)*" 1 , (18) 

As was explained above, the matrix \Pi is a specification 
of ^ corresponding to a particular value of the matrix 
spectral parameter: 

'AS 
£i = | AS 
A 

We denote as $o the fundamental matrix of solutions 
for the linear system Eqs.Q,© f° r A = Ao. It can 
be shown that for the value of the spectral parameter 
A = \q the fundamental matrix is <J>o = (^o 1 V ■ Since 
the subspace of solutions corresponding to Aq is two di- 
mensional, the matrix *i is constructed as follows. The 

(3) (1) (2) (3) 

vector = ci$Q + c 2^o + c 3®o ^ s a general so- 
lution of the linear problem with A = Ao. Here upper 
index in the brackets i — 1,2,3 denotes a vector-column. 
To satisfy the structure of the operator 5 Eq. fTTjl we re- 
quire that (^f\ ' , ' ') = 0, where (•, ■) denotes a scalar 
product of two vectors in 3D complex vector space, and 

(1 2) 

the vectors ' correspond to A = Aq. Due to the 
definition =6ij. 

We can easily find two appropriate orthogonal vectors 



(1,2). 



(1) 



(4 + cl)^ 1} -cl(^>+^') 



(2) 



s(3h 



(2) _ 



4% 



(2) 



(3) 



The algorithm for finding new solutions of the nonlinear 
system Eqs.0,® can be formulated as follows: Find a 
solution $ °f the associated linear system Eqs.0,0, 
corresponding to a certain " seed" solution of the nonlin- 
ear system Eqs.©,©; Build and build 5(A), then 



use the dressing transformation Eq. (|16fl . It is straightfor- 
ward to show that for the state Eq.© of the atom- field 
system a general solution of linear system Eqs.Q,© can 
be represented in the following form 



$o = 



T(t, A) 



(19) 



W(r,A) 



where a 2 x 2 matrix T(r, A) is defined through two com- 
plex functions w(t, A) and z(t, A) as follows 

T(t, A) = (I + W(t, A)) e z( - T ' X \ (20) 

-w*(t,X) 
w{t, A) 

\ -i|t + z*(t,A) 

Here I is a 2 x 2 identity matrix. The function w(t, A) 
satisfies the Riccati equation 

-id T w(r, A) = -Xw(t, A) + ifi( T ) - ( r ) w 2 (t, A), 

(21) 

and the function z(t, A) is defined through w(t, A): 

-id T z(T,\) = ~Sl*(T)w(T,\). 



(22) 



It is easy to check that $o has the same form as $o with 
A replaced by A* , and T replaced by 



T(r,A) 



(1 + W(r, A*))e- 



-Z*(t,A*) 



1 + w(t, \*)w*(t, A*) 



(23) 



Applying the procedure described above, we find for the 
fields 

n« = -2 3(0)3,1 = (24) 
-2(A - \l)e- l ^ {w ^ x) % 2+eV3 \ 

fj 6 =fi(r)- 2 5(0)3,2 = (25) 

n(r) - 2(A - K) ^-MrX)e^r,X)e-2 + e^) 

The corresponding density matrix p = \tp a t)(ipat\ reads 

1,7, \ _ 5(A) 1 , 1 |l)+H(A) 2 ,i|2)+H(A) 3 ,i|3) _ (9r ^ 
Wat) ~ [a-TaI - \ M ) 

( _j_ (Ao-Ag) \ m\ , 

V|A -A| ^ \\o-A\AT J I 1 / ^ 

e~' V1 (Ap-Ag) ( (e^2_^«( T; A) e ^3) (w^A)e^+e^3) , . \ 
[Ao^Aj ^ 7\7 I 2 ? + A? \ S I ) ■ 

Here, ipi is the phase of the coefficient c\. The module 
of this coefficient can be set to unity without loss of gen- 
erality. For shorter notations we have also defined two 
phases (cf. Eas-fEty. JU): 



vg C 
A-A 



(p 2 = Z(r, A)i 4 + log(c 2 ) - | (\t 

= z ( T ,A)+log( C2 )- 5 ^ y , 
¥>3 = Z(t, A) 2j2 + log(c 3 ) - | (At + ^ 
= -zAr + z*(r, A) +log( C3 )-^ y , 



(27) 
(28) 
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and the normalization function 

AT = 1 + Re [(w(t, A) - w(t, X*))e^+^] 

+ (1 + \w(t, X)\ 2 )(e V2+ ^ + ef 3+ ^). (29) 

In conclusion of this section, we note an important differ- 
ence between y>2 and ip^. It will be shown below that for a 
constant or slowly varying background field the function 
z(t, A) is of the same order of magnitude as the control 
field intensity |i7(r)| 2 . Therefore, for small intensities the 
phase if2 is slowly varying in time and describes the slow- 
light soliton, while the phase <£>3 is varying with a speed 
close to the speed of light in the vacuum due to the term 
At. 



III. 



THE TRANSPARENCY GATE 



In this section we introduce a concept of fast and slow- 
light solitons in the A-mcdium and explain how the non- 
linear interplay between the solitons leads to a possibility 
to control transparency of the medium. We discuss first 
the mechanism of transparency control for the slow-light 
soliton. We explain how the fast soliton propagating in 
the a channel hops to the b channel where the slow-light 
soliton is propagating. The fast soliton then destroys 
the slow-light soliton, thus stopping the propagation of 
the latter, and then disappears itself due to the strong 
relaxation in the system. 

As was indicated above, in this work we consider exact so- 
lutions of the Maxwell-Bloch system Eqs.©,© existing 
on some finite background. The background field plays 
the same role as the controlling field in the conventional 
linear theory of EIT, but it enters the exact solutions 
as a parameter in a substantially nonlinear fashion. We 
start with the case of a time independent field specified 
as follows 



(30) 



Here k <C k a ^ is introduced in order to take into account 
small spatial variations of the phase. The intensity of 
the background field Slo is an experimentally adjustable 
parameter, which provides control over the transparency 
of optical gates and determines the speed of the slow- 
light soliton. The Maxwell-Bloch system Eqs.©,© is 
satisfied with the following initial state of atoms 



1 




i^oe^ £(§-A) 



(31) 



The parameter x determines the population of the ex- 
cited state and has to be larger than 2A. It is important 
to notice that for a time-independent background field 
atoms can be prepared in a mixture of dark state and po- 
larized states only for nonvanishing parameter k, which 
allows to access a wider range of physically interesting 
situations. 



For a time-independent background field we immediately 
find solutions of Eas. (|2~l")) . (|2*2l . viz. 



w(t, A) = w 



z(t, A) = z t ee |O e lkt: w T 



i £7 n r 



2(A+VA 2 +O 2 ) 



(32) 
(33) 



As we noted in the previous section, when fl b depends 
on ( 



(*o) 



i,i 



= A Xr+ 



A- A . 



and the structure of the solution T(t, A) in Ea. l|2U[l is 
slightly modified, i.e. we have to replace Z(t, A) with 



Zi(r, A) = Z(t + 



where 173 is the Pauli matrix. Hence, the phases Eas. (|27ll . 
Pfy read 



ip 2 = log(c 2 ) 



i n t i k £ 

2(A+ V /A^+SV) 2 



i(3kx+2ky/ A 2 +Q 2 -2i/ )C 
4(A-A) - 



1 f \ i\r i\/A 2 +f2o 2 T ikC 
Lp 3 = l0g(c 3 ) - — Y g + 2 

i(3kx-2ky/ A 2 +S1 2 -2v K 
4(A-A) • 



(35) 



(36) 



Using the general solution Eas.Q24 [l .H25 |) we can find the 
dynamics describing the formation of the transparency 
gate for initial conditions specified in Eas. (|30|) . H31|l . For 
simplicity, in this section we take the spectral parameter 
to be purely imaginary, Ao = i eo, and for a solitonic type 
of solution eo > ^o- The solution corresponding to the 
phases Eas. (|35|l . I|35|l describes the nonlinear interaction 
of fast and slow- light solitons. This solution is parameter- 
ized by the constants C2.3 defining the position and phase 
of the two solitons. As we have already indicated above, 
the phase Lp 2 determines the position of the slow-light 
soliton whereas y? 3 determines the position of the fast 
signal. In practice these constants 02,3 are defined by the 
initial condition, which specifies the actual pulse of light 
entering the medium at the point z = 0. To understand 
the structure of the slow-light soliton one can set C3 = 0. 
This choice corresponds to taking the fast soliton to -co 
in the variable r. Indeed, this specification removes the 
fast pulse component corresponding to ^3 and thus sin- 
gles out the slow- light soliton part. The slow- light soliton 
solution assumes then the following form: 



(Ag-A) w e'' Im w-w» 



where 



&b = -£loe ik< - tanh^) 



Re^ + -log(l + |w | 2 ) 



(37) 
(38) 

(39) 
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is the phase of the slow-light soliton. For simplicity, in 
the following we let k = 0. From the expression above 

and in the simplifying approximation « 1, A = 

£ 

the group velocity of the slow-light soliton can be easily 
derived: 



' v ' 



(40) 



The pure state of the atomic subsystem corresponding to 
the slow-light soliton solution reads 



\i>at 



RcA— A— ilmA tanh t 



|A-A| 



ID 



S2„ 



2|A-A|u> 



|2) 



"2|A-A| 



13). 



(41) 



Notice that the population of the upper level |3) is pro- 
portional to the intensity of the background field. The 
speed of the slow- light soliton is also proportional to Slg. 
This means that the slower the soliton, the smaller the 
population of the level |3) and, therefore, the dynamics 
of the nonlinear system as a whole is less affected by the 
relaxation process. 

To understand the structure of the fast soliton one can 
choose C2 = 0. We then arrive at an expression describing 
a signal moving on the constant background with the 
speed of light (fast soliton): 



A (A*-A)e i ( Imv 3- ,p i ) , , , , 

n b = -Q e lk ( tanh(0/), 
where the phase of the fast soliton is 

<j>f = Reip 3 + -log(l + K| 2 ). 
The atomic state is described by the function 



(42) 



Hat) = 



RcA — A — ilmA tanh * 



|A-A| 



e%i2) 



2|A-A| 



2|A-A| 



|3>. 



(43) 



We emphasize the principal difference between fast and 
slow-light solitons. The slow-light soliton vanishes when 
the controlling field is zero due to the factor wo in 
Ea. (|37fl . Another important feature is that for the slow- 
light soliton the population of the upper level |3) is pro- 
portional to Oo and is small for small background fields 
and therefore stable to optical relaxation. In contrast, 
the amplitude of the fast signal is not limited by Qq and 
is determined by the spectral parameter sq. The popu- 
lation of the level |3) is also defined by sq, which means 
that for large spectral parameters the fast signal will be 
attenuated by the relaxation. As we discuss below, in 
the absence of the background field the fast signal be- 
haves as a conventional SIT soliton in a two level system 
|1><— > |3>. 

Figure [21 illustrates the propagation and collision of 
the fast and slow-light solitons according to equation 



Eas. l|24(l . H25|) . The figure for I a shows the intensities 
of the signals in channel a. We see that before the colli- 
sion only the slow-light soliton exists in channel a, while 
after the collision the slow-light soliton disappears and 
a fast intensive signal appears, whose velocity is slightly 
below the speed of light. The figure for I a is comple- 
mented by the figure for the intensity If, of the field in 
channel b. The slow-light soliton corresponds to a groove 
in the background field £Iq. It is clearly seen that after 
the collision the slow-light soliton ceases propagating in 
channel 6, while some trace of the fast soliton still can 
be noticed in that channel. The process described above 
can be summarized as if the fast soliton destroys the slow- 
light soliton. The notion of a transparency gate requires 
the existence of two distinctly different regimes, which 
are transparent (open gate), and opaque (closed gate). 
In the absence of the fast soliton the gate is open for the 
slow-light soliton. When the fast soliton is present the 
slow-light soliton is destroyed, while the fast intensive 
signal created after the collision in channel a is attenu- 
ated due to strong relaxation in the atomic subsystem. 
The gate thus closes in the course of the dynamics due 
to the relaxation process. To further explain this process 
we provide the Fig. [2] plots for populations of the levels 
1 2) and |3>. 

Notice that before the collision the population of the up- 
per atom level |3) is negligible and is approximately given 
by the formula for the slow-light soliton solution Eq. I|41|) 
(see the lower right plot of P3). The populations of the 
lower levels |1,2) are determined by the slow-light soli- 
ton (see the lower left plot of PI). Indeed, the fast signal 
existing in channel b does not interact with the atoms 
because at the onset of the dynamics their state coin- 
cides with the dark state |1). Figure |3 shows that after 
the collision the atoms of the active medium are highly 
excited and therefore the level |3) is strongly populated. 
This leads to the fast attenuation of the rapid intensive 
signal in channel a due to the relaxation. The optical 
gate closes. 

To this point we have described a mechanism of con- 
trolling the transparency of the medium for a particular 
type of slowly moving signals. We now discuss a possi- 
bility to read information stored in the atomic subsys- 
tem. Let us assume that the background field vanishes, 
i.e. Slo = 0. As was explained above, the speed of the 
slow-light soliton then vanishes as well. However, the in- 
formation about polarization of the slow-light signal is 
stored in the atomic subsystem. This effect can be in- 
terpreted in terms of the concept of a polariton, which 
is a collective excitation of the overall atom-field system. 
The notion of a polariton for the A-system has been used 
before. In the linear case the dark-state polariton was 
discussed in |3fj |. In the strongly nonlinear regime, which 
is the case for the present work, a similar interpretation 
is possible. Indeed, the field component of a slow-light 
soliton solution can be interpreted as the light contribu- 
tion into the slow-light polariton. When the controlling 
field Qq vanishes this contribution also vanishes, along 
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FIG. 2: Knocking down the slow-light soliton. The two upper 
plots correspond to the dynamics of the intensities of the fields 
fi a and Qb- The two lower plots show the populations of the 
levels 1 2} and |3). The parameters of plotted solutions are 
C2 = C3 = 1, Ao = 4.1i, A = 0. 

with the speed of the polariton. The latter then contains 
only excitations in the atomic subsystem. The general 
solution Eas. l24|) . I|25[) is then reduced to the form: 

_ _ 2ie exp[i 2( ^ 2) +ilog( 1 ^ T )-»y 1 | 
a ~ cosh(0 sO ) + | cxp[20 /o -0so] ' 

n b = e llpl+log(c2) -^T^ fi a , (44) 

where 4> s q — ,, + log( Ic2 1 ) is the phase of the slow- 

light soliton, and <f> f0 = e q t + 2( ^'^ 2) + log(|c 3 |) is the 
phase of the fast soliton for the vanishing background Slo- 
The form of the fields resembles a superposition of fast 
and slow-light solitons in Eqs.JS^J, (|25|l with the vanish- 
ing velocity of the slow-light soliton. This last exponen- 
tial term in the denominator represents the fast signal 
contribution. The component containing hyperbolic co- 
sine describes the information about the slow-light soli- 
ton stored in the medium after the soliton was stopped. 
The overall picture of dynamics corresponds to the scat- 
tering of the fast soliton on this localized atomic polariza- 
tion pattern. The atomic state describing this scattering 
reads 

w = -( v^S:^ + 7^)' 1 > (45) 

For C3 = the fields vanish, while the atomic state re- 
duces to a form corresponding to a stopped polariton 
described by Ea. (|41|l with Q = 0. In other words, when 
the slow-light soliton is completely stopped, the informa- 
tion borne by the soliton is stored in the spin polarization 




FIG. 3: Reading the optical information by the fast soliton. 
The two upper plots illustrate the dynamics of the fields Q a 
and Qt- The two lower plots show the populations of the 
levels 1 2} and j3). The standing peak on the plot for Pi corre- 
sponds to the stored information in the form of the localized 
polarization. The rapidly moving localized excitation of the 
atoms given on the plot for P3 represents the act of reading. 
The background field Qo = 0, and the coupling constant is 
the same as before, i.e. z^o = 4.5. 

of the atoms. As long as the upper state |3) is not popu- 
lated, the state of the atomic subsystem is not sensitive 
to the destructive influence of the optical relaxation pro- 
cesses. 

The conventional way 0] to read the information stored 
in the atoms is to increase the intensity of the background 
field. Our method of reading the information is different. 
We propose to send the fast soliton into the space domain 
in the active medium, where the information is stored. 
The polarization in the domain is then flipped by the fast 
signal. This is how the reading of information is realized. 
This way of reading optical information is advantageous 
because it involves fast easily detectable processes. Fig- 
ure |3| illustrates the mechanism of the reading. Notice 
that the act of reading, based on the polarization flip- 
ping induced by the fast signal, can be realized on a very 
short time scale compared to typical relaxation times. 

IV. NONADIABATIC MANIPULATION OF 
SLOW-LIGHT SOLITONS 

In this section we discuss an exactly solvable, though 
physical, case describing controlled preparation, manip- 
ulation and readout of slow-light solitons in atomic va- 
pors and Bose-Einstein condensates. The group velocity 
of the slow-light soliton depends explicitly on the field 

Qo i- e - v g w c^- (cf. Eq.JJJJ). This expression im- 
mediately suggests a plausible conjecture that when the 
controlling field is switched off the soliton stops propagat- 
ing while the information borne by the soliton remains in 
the medium in the form of a spatially localized polariza- 
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tion pattern, i.e. optical memory, which can be recovered 
later. For brevity we refer to this pattern as to "memory 
bit". 
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FIG. 4: The intensity of the field Qb as the function of time 
t. The thick solid curve shows the time-dependence of the 
controlling background field at entrance into the medium, i.e. 
at z = 0. We also plot \Clb\ 2 at z = 6 (dashed curve) and 
at z = 12 (dotted curve). We choose a = 4, and the delay 
interval is T — T\ — 3. The dimensionless units are defined in 
the text. 

We consider the following scenario for the dynamics (see 
FigQJ. Before r = we create in the medium a slow- 
light soliton, and assume it is propagating on the con- 
stant background ilo- We then slow-down the soliton 
by switching off the laser source of the background field. 
We assume an exponential decay of the background field 
with a decay constant a, i.e. Qq£ ~ aT . At a certain mo- 
ment of time, say, Ti = 4 /a the field becomes negligible. 
Therefore, we cut off the exponential tail and approxi- 
mate it by zero. At this step the soliton is completely 
stopped. The position where the soliton stops, depends 
on the decay constant a and on the moment when we 
switch the laser off. The information borne by the soli- 
ton is stored in the form of spatially localized polariza- 
tion. This formation can live a relatively long time in 
atomic vapors or BEC 1 1 j . At the moment T we re- 
store the slow-light soliton by abruptly switching on the 
laser. The whole dynamics is divided into four time in- 
tervals ULo^ = (-oo,0]U[0,Ti]U(Ti,T]U(T,oo). The 
time-dependence of the intensity of the background field 
at entrance into the medium is given in Fig^J 
Before the soliton enters the medium, the physical sys- 
tem is assumed to be prepared in the state described by 
Eq.(PJ). The function £1(t) now models switching the con- 
trolling field off and on again. This function reads (cf. 
Fig0: 

fl(r) = n [0(-r) + e-^(6(r) - 6(r - Ti)) (46) 
+6(r-T)]. 

Here O(-) is the Heaviside step function with 9(0) = \. 
For the state Eq.@ with Ea. H46|) we exactly solve the 
nonlinear system Eqs.JSJ),!© as well as the auxiliary scat- 
tering problem Eqs.Q, (JSJ) underlying its complete inte- 
grability. The latter result is the cornerstone of analyti- 
cal progress achieved in this section. This result allows to 



mounting a soliton on the background Eas.(jfj|l. (|46(l using 
the Darboux-Backlund transformation described in sec- 
tion^ (cf . 24] ) . According to the results of sectionlTTlthc 
one-soliton solution corresponding to the time dependent 
background Eq.JSJ) reads 



n h = 



V1 + Ht,A)| 2 
_ (A-A>(t,A) j, 



(47) 



e' ps sech0 s - O(r), 



1 + Ht,A)| 2 
with the atomic state p = \ipat)(tpat\, where 

I 7 \ RcA— A— ilmA tanh <f> s i -i \ 

Wat) — | A- A | I / 

+ 2|A-AMt,A)I 2 ) ~ 2 1 A | I 3 )' 



(48) 



Here. 



o + ^lm^^+Re(z(r, A)) 



In Vl + Kt, A) I 
2 



n - ^-Re - 1 - + Im(z(r, A)), 



A- A 

where A is an arbitrary complex parameter. The func- 
tions w(t, A), z(t, A) are of piecewise form, specific to 
each time region T>i. For clarity we organize elements 
of the solution corresponding to different time regions 
in Table [I] The function wo is defined as in Eq. IMl'l) 
with k = 0, the index of Bessel functions is defined as 
7 = (a + i\)/(2a). The values Ci of the constant C for 
each time region T>i are specified in the rightmost col- 
umn of the table, the moment of time T is chosen as 
in FigH i.e. T = A/a + 3. Notice that in the table 
U>2 = Wi(oc, A) and z-i = Zi(oo, A). Therefore the solu- 
tion in the region T>2 is parameterized by the asymptotic 
values of the data for the region T>\ corresponding to the 
absence of cut-off of the exponentially vanishing tail. The 
region T>2 describes the phase when the slow-light soliton 
is stopped: the fields vanish, while the information borne 
by the soliton is stored in the medium in the form of spa- 
tially localized polarization. At the time T the laser is 
instantly turned on again. The stored localized polar- 
ization then generates a moving slow-light soliton. This 
process is described by the solution in the region T>3 . Ex- 
cept for the moment Ti, the functions w, z are continuous 
in t. This ensures that the physical variables such as the 
wave-function and field amplitudes evolve continuously. 
We demonstrate typical dynamics of the intensity of the 
fields in Figs^and|S| In FigQJthe decaying shock wave, 
whose front has an exponential profile, propagates with 
the speed of light, reaches the slow-light soliton and stops 
it. An intense and narrow peak is developing in the back- 
ground field when the shock wave hits the soliton (dotted 
curve). This peak signifies a transfer of energy from the 
soliton to the background field. After the auxiliary laser 
has been switched on again, another step- like shock wave 
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TABLE I: Exact analytical solution 



V 


0(t) 


u>(r, A) 




c 

















. CJ 1— r( - Mr^) _J 7-l( _ TS~") 




-tooJ7(-^)+^7-l(-^) 






J 1 _ 7 (-^) +iTO0 J_ 7 (-^) 


v 2 








, c (-^) 7 /r(i-7) 


C 2 =Ci 


£>3 


n 


n tan(i^A2+n2 (T _ T) ^ 


j Ce 2 +e 2 




A tan i ^/A 2 +D 2 , (t-T) J -i,/A 2 +n;i 


111 C+l 

+Z2 
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FIG. 5: Contour plot of the intensity of Cl a . We choose 
A = — 4.1i and zero detuning, A = 0. The break-up area in 
between the two solitonic trails manifests the creation of a 
standing memory bit in the medium. 



reaches the localized polarization formed in the medium 
by the incoming soliton, and retrieves stored informa- 
tion in the form of a new slow-light soliton. A narrow 
and deep depression in the intensity plot means now the 
energy transfers in the opposite direction, i.e. from the 
background field to the restored soliton. The dynamics 
of the field £l a is plotted in Fig|5J The contour plot shows 
that in the process of rapid deceleration the solitonic trail 
profiles end sharply. Notice that the characteristics of the 
restored pulse, i.e. the width and group velocity, are very 
close to those of the input signal existing in the medium 
before the stopping. 

We now calculate the half-width of the polarization flip 
written into the medium after the soliton is completely 
stopped. It reads 



W s =4cln(2 + \/3) 



|A~A| 2 
u \Im(X)\ 



(49) 



It is important to notice that the width Ea. H49|) of the 
optical memory formation does not depend on the rate a. 
In other words, the width of the memory bit is not sensi- 



tive to how rapidly, i.e. nonadiabatically, the controlling 
field is switched off. This leads to an important conclu- 
sion. Indeed, through the variation of the experimentally 
adjustable parameter a, it is possible to control the lo- 
cation of the memory bit, while its characteristic size re- 
mains intact. Dutton and Hau have already reported |Tl| 
that when the switching is made quickly compared to the 
natural lifetime of the upper level, the adiabatic assump- 
tions are no longer valid, but, remarkably, the quality of 
the storage is not reduced in the nonadiabatic regime. 
Our analytical result is in excellent agreement with this 
observation. 

The group velocity of the slow-light soliton reads 



c 



Mt,a)| 2 



Mi+Mt,a)|2) 

2|A-A|2 



\w(t,XW 



(50) 



Notice that in the case of the constant background field, 
i.e. in the case a = 0, the conventional expressions for the 
slow-light soliton Eas. (|37(l . l38|l along with the expression 
for the group velocity Ea. (|40|) - the main motivational 
quantity for this report - can be readily recovered from 
Eqs.j4 T) l.(l5C) |) . 

The distance C s (a) that the slow-light soliton will propa- 
gate from the moment when the laser is switched off until 
the full stop, is 



C s {a) 



2c|A-A| 2 
t/ Im(A) 



I T — OO 

It=0 



2e|A-A| 2 
"o|Im(A)| 



lnyi + K| 2 -Rc(z(oo,A)) . (51) 



It is evident from our solution Eq.JSlJ that the soliton 
possesses some inertia or, in other words, a momentum of 
motion. Indeed, even if the controlling field is switched 
off instantly (notice that lim Rez(oo, A) =0), the soli- 

a — >oo 

ton will still propagate over some finite distance after the 
shock wave of the vanishing field, propagating with the 
speed of light, has reached the soliton. 
In FigElwe show the dynamics of localized polarization 
corresponding to soliton dynamics shown in FigEl The 
plot is remarkably similar to the corresponding figure in 
reference describing recent experiments. The central 
part of the plot in FigEl shows the standing memory bit 
imprinted by the slow-light soliton. Notice that in the 
presence of the soliton the population flip from level |1) 
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FIG. 6: Population of level |2). Here, A = -4.1i and A = 
0. The time interval 1 < t < 4 corresponds to a standing 
localized polarization flip. Compare with experimental results 
in Fig.l of ref. Q. 



to 1 2) in the center of the peak is almost complete. This 
property of the solution manifests the major distinction 
between the strongly nonlinear regime considered in our 
paper and the linear EIT theory. As was already pointed 
out earlier 0, 0, |2(j, the regime of two fields being 
comparable in magnitude opens up new avenues for an 
effective control over superposition of two lower states of 
the atoms. Changing the parameters of the slow-light 
soliton we can coherently drive the system to access any 
point on the Bloch sphere, which describes the lower lev- 
els. For zero detuning, the solution discussed here shows 
that when the field vanishes the maximum population of 
the second level reaches unity. Using Eq. (|4*S|) it is also not 
difficult to estimate the maximum population of the level 
1 2) for finite detuning after the soliton was completely 
stopped: |A| 2 /|A — A| 2 . Notice that only a small fraction 
of the total population is located in the upper level |3) 
and provides for some atom-field interaction. This popu- 
lation is proportional to |f2o| 2 - Numerical studies of the 
Maxwell-Bloch equations with relaxation terms included 
[TlllT^ I show that for experimentally feasible group veloc- 
ities of the slow-light pulses, i.e. when the maximum in- 
tensity of the controlling field is not very high, the pulses 
are stable against relaxation from level |3). Here we con- 
sider the same range of parameters. Therefore, the de- 
structive influence of relaxation on our solution Eq. 14711 
is negligible. 



V. THE CASE OF ARBITRARY 
CONTROLLING FIELD 

In this section we build a single-soliton solution on the 
background of the state of the overall atom-field system 
described by Eq. Q for a quit general (complex) control- 
ling field The single-soliton solution corresponding 
to the background field fi(r) is given by Eas. p7Jl . (|4"g|) 



with the functions w(t, A) and z(t, A) defined below. 
We envisage the dynamics scenario similar to that of the 
previous section. We assume that the slow-light soliton 
is propagating in nonlinear superposition with the back- 
ground field, which is constant at r — > -co and vanishes 
at r — > +co. The speed of the slow-light soliton is con- 
trolled by the intensity of the background field. There- 
fore, when the background field decreases, the slow-light 
soliton slows down and stops, eventually disappearing 
and leaving behind a standing localized polarization flip, 
i.e. optical memory bit. Should the background field 
increase, the soliton will emerge again and accelerate ac- 
cordingly. 

To be specific, we define the asymptotic behavior for the 
field f2(r) in the form 



flOr 



-co) = f2oj — ► +co) = 0. 



(52) 



The asymptotic boundary conditions Eq. I|52|) dictate the 
following asymptotic behavior for the functions w(t, A) 
and z(t, A) defined by Eqs. l|2T J) .ip2 )l : 



w(— oo, A) = wq = 
z{— co, A) 



Z T 



2k(X) 
.Ifiol 



4fc(A) ' 



, w(+co,A) = 0, (53) 
(54) 



where k(X) = (A + ^JX 2 + \n \ 2 )/2. The function z(r, A) 
satisfying the asymptotical conditions Eq. (|54|) reads 



z(t,X) = z q t + I [-n*(T')w(T',X)-Z )Q(T-T')dT' 



The function w{r, A) is defined by the relations 



(55) 



w(t,X) = i J e~ lk{T - s) ®{ T - s)w(s,X)ds, (56) 

— OO 

CI(t) 1 



w(t,X) 



Here 0(r) is the Heaviside step function. We rewrite the 
relations Eas. l|56|l . (|57(l in the form of nonlinear integral 
equation, viz. 

w(t,A) = J e- tk ^-^e{T-s)w{s,X)ds 

— OO 

OO / * 2 \ 

. J e -ifc(T- s ) ( T - s ) rairltf^A)- l^Lj ds. (58) 

Hence, we can construct a solution w(r, A) iterating 
Eci. (|55|) and starting iterations from w (r, A) = |^(r). 
Notice that the last term in Eq. I|57ll provides a correction 
of order k~ 2 , because the function w(t, A) asymptotically 
behaves as 1/fc. In the adiabatic regime when the back- 
ground field varies slowly, i.e. all derivatives of ^(r) are 



11 



much smaller than k, we can integrate Ea. (|56|l by parts. 
Preserving only the lowest order term with respect to k 
we obtain 



w(t, A) 



2k ' 



(59) 



Hence, at the lowest order in k we find 



z(t,A)«z t+ / [ ^\Q(t')\ 2 - z ) e(T-r')dr' . (60) 



As one can easily observe this expression is in agreement 
with the asymptotic condition Ea. (|54() . 
For an arbitrary dependence of the background field on 
the retarded time r, the speed of the slow-light soliton 
can be represented in the following form: 



d T > 



c d T <f> s — i 

It can be readily seen that 

d<f> s = Im(A)| W (r,A)| 2 = 
dr l+\wW,X)\ 2 ' ac 



v ° T 

y Im r 



(61) 



(62) 



We have thus found a general solution for the velocity 
v g of the slow-light soliton propagating on an arbitrary 
time-dependent background field in terms of the function 
w(t, A) given by Eq.JSHJ- This result provides a new way 
to study dynamics of localized optical signals in the non- 
linear EIT systems. It allows us to easily suggest different 
schemes to slow-light down, stop, and reaccelerate slow- 
light solitonic contribution in the probing pulse. With 
such techniques one can introduce a concept of probing 
different regions of the media by changing the time that 
the soliton dwells at a particular location. This time 



is important in situations where the interaction between 
light and some impurities inside the EIT medium is weak 
and requires slowing the signal down in the vicinity of 
these impurities in order to gain more information about 
the structure of the medium. 

We also introduce a notion of the distance £[f2] that the 
slow-light soliton will propagate until it fully stops. This 
quantity is important because it describes the location 
of an imprinted memory bit. The brackets [•] indicate a 
functional dependence of the distance on the controlling 
field f2(r). To begin with we consider the case when the 
field is instantly switched off at the moment r = 0, i.e. 
£1(t) = ^qQ(-t). Then we easily find the solution for w 
and z: 



w(t, A) = w (6(-r) + 9(r) e~ aT ) , z(r, A) = z @(-t)t. 

Hence, we can obtain the distance Cq that the soliton 
will propagate from the moment r = until its complete 
stop at r — > oo: 



Co 



c|A- A| 2 



ln(l 



\ w o\ 



u \lm(X)\ 

Here we make use of the assumption that Im(A) < 0. 
Now, we can give the definition of the distance £[f2] for a 
generic field f2(r) satisfying the conditions Eq.(j52J|. It is 
convenient to define it as a relative distance, namely the 
difference between the absolute coordinate of the stopped 
signal at the maximum of the signal and the distance Cq. 
The relative distance reads: 



C[Cl] = 



2e|A - A| s 
u lm(X) 



Re 



-Q*(r)w(r,X)- z e(-T) ) dr. 



Using the representation Ea. H57|) we find 



-f-OG -j-OO 

£[»] = 2c| t~ A f Re / f e - ,fc(T - s) e(r- S )fMeH^»( S ,A)]^rfr| (63) 



f Im(A) 



If we assume that O(r) is a smooth function and sub- 
stitute the solution for w(t, A), we find the result in the 
form of a series 

2c|A-A|% f^I n \ 

where are regularized Zakharov-Shabat 

functionals |2j. The first two functionals 

read h[Sl] = - (|tt(r)| 2 - |tt | 2 e(-r)) dr, 

W = mI-oo( n *( s ) 9 Ms) - tt(s)d s n*(s))ds. The 
other functionals can be obtained through the iter- 



ation procedure described above. As it is usual for 
the boundary conditions of finite density type, I\ is 
not a proper functional on the complex manifold of 
physical observables, in the sense described in |21|. In 
that sense all other functionals in the expansion with 
respect to k are proper. It is a plausible conjecture that 
the minimum of the functional of length, Ea. l63|) . i.e. 
5jC[Q]/SQ = with S 2 C[n]/Sn 2 > 0, is achieved when 
the controlling field is switched off instantly. Therefore it 
seems intuitively correct that the minimum is delivered 
by the the function fioO(— r) discussed above. This 
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conjecture is also supported by the case discussed in 
section HVI when the controlling field vanishes exponen- 
tially, i.e. Q(t) = n (6(-r) + 6(r)e- Qr ). In this case 
the minimum of length is delivered by a singular limit 
a — ► oo, i.e. in the regime of instant switching off of the 
controlling field. 

Another important characteristics of the system is the 
shape of the imprinted signal. It is easy to show that the 
width Wo of the imprinted memory bit is not sensitive 
to the functional form of f2(r) and is the same as given 
by Ea. (|49|l . In other words, this exact result is valid 
regardless of how rapidly we switch the background field 
off. This means that specification of f2(r) only influences 
the location of the stored signal and not its shape. This 
result is strongly supported by recent experiments |llj| . 
This reference emphasizes the phenomenological fact that 
the quality of the storage is not sensitive to the regime 
of switching off of the control laser. Our exact result 
Ea. (149(1 provides a rigorous basis for this experimental 
observation. 

To conclude this section we discuss the applicability of 
the concept of effective time to the regime of nonadia- 
batic variations of the controlling field. This concept was 
used before in [il lhillicj for approximative descriptions of 
pulse propagation on the background of a time-dependent 
controlling field. To account for this dependence, these 
references introduce an effective time variable (a "scaled 
time" of ref.B and see the function Z{r) of ref. 14] ). 
The reference jlj|) shows the concept of effective (scaled) 
time to be very useful in the regime of linear EIT, while 
ref. Q demonstrates applicability of this concept in the 
strongly nonlinear, though adiabatic, regime. The ef- 
fective time approximation, in our formulation, is given 
by expressions Eas.ljK ^ .lpI Hl employed in the slow- light 
soliton solution Eas. (|17|l . l|4^|l . In Fig. [7] we compare the 
resulting approximate solution against an exactly solv- 
able case discussed in section IIVI We point out that 
the method of effective time has a rather limited ap- 
plicability in the strongly nonlinear noadiabatic regime, 
which is most interesting for modern experiments |ll| . 
Indeed, Fig. demonstrates that this method as applied 
to the slow-light soliton gives rise to large errors in the 
regime of nonadiabatic dynamics. Therefore heuristic at- 
tempts to substitute an effective time into the slow- 
light soliton solution Eas. (|37f) . (|^8 |l is largely not accurate. 



VI. CONCLUSIONS AND DISCUSSION 

In this paper we have discussed a physically realistic, 
exactly solvable model of manipulation, i.e. prepara- 
tion, control and readout, of optical memory bits in the 
strongly nonlinear, and more importantly, nonadiabatic 
regime. We discussed a concept of the transparency gate 
for the slow-light solitons in a A-type medium. We ex- 
plained how the fast soliton can destroy the slow-light 
soliton and close the gate for the latter. We also de- 
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FIG. 7: Here e = 4.1, A = 0, and fi(r) = 0.5e" 4T . It is 
clear that when the background field changes relatively fast 
the exact solution given by solid and dotted lines are differ- 
ent from approximative one (dashed and dot-and-dash lines, 
correspondingly) by up to 40%. 



scribed the process of reading optical information, writ- 
ten into the active medium by the slow- light solitons. 
We have investigated a mechanism of dynamical control 
of the slow-light soliton whose group velocity explicitly 
depends on the background field. For a quite general 
background field, we found the location and shape of the 
memory bit written into the medium upon stopping the 
signal. Remarkably, the width of this spatially localized 
standing polarization flip is not sensitive at all to the 
functional form of the controlling field and is defined by 
the parameters of the slow-light soliton only. 
It is worth discussing here a possibility to actually cre- 
ate in the A-type atomic medium the slow-light solitons. 
The general physical feature underlying the mathemat- 
ical property of complete intcgrability is a delicate bal- 
ance between the effects of dispersion and nonlinearity 
inherent in the medium. Provided that this balance is 
observed and the system is completely integrable, it is 
a general fact that virtually any sufficiently intensive 
localized initial condition creates solitons. The overall 
picture of nonlinear dynamics can be roughly described 
as follows. The evolving signal created by the incident 
pulse in the course of the dynamics separates into a num- 
ber of solitons and a decaying tail. The latter vanishes 
in due course. The soliton-like signals survive (ideally, 
i.e. in the absence of dissipation) for infinitely long time. 
When the physical conditions underlying the complete 
integrability of the optical system are met, the general 
picture of the nonlinear dynamics is similar to that de- 
scribed above. Namely, a fairly arbitrary localized and 
intensive initial signal creates in the course of nonlinear 
dynamics a number of slow-light solitons. The number of 
slow-light solitons can be derived from the analysis of the 
corresponding zero-curvature representation. The signal 
Eas. l|37|l . (|38|) is a generic soliton-like solution of the non- 
linear system Eqs.lJSJ), © and therefore it is very plausi- 
ble that such signal can be created. To conclude, we want 
to emphasize that in our considerations the distinguished 
role is assigned to the background field f2( T ) that turns 
out to be a nonlinear analog of the conventional control- 
ling field appearing in the linear EIT formulation. The 
difference between the linear and nonlinear cases lies in 
the fact that in the nonlinear case the control field and 
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the slow-light soliton solution are present in the same Work at Los Alamos National Laboratory is supported 
channel in an inseparable nonlinear superposition. by the USDoE. 
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